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ABSTRACT 



Aims. We aim at getting high spatial resolution information on the dusty core of bipolar planetary nebulae to directly constrain the 
shaping process. 

Methods. We present observations of the dusty core of the extreme bipolar planetary nebula Menzel3 (Mz3, Hen 2-154, the Ant) 
taken with the mid-infrared interferometer MIDI/VLTI and the adaptive optics NACOA'LT. 

Results. The core of Mz 3 is clearly resolved with MIDI in the interferometric mode, whereas it is unresolved from the Ks to the N 
bands with single dish 8.2m observations on a scale ranging from 60 to 250mas. A striking dependence of the dust core size with 
the PA angle of the baselines is observed, that is highly suggestive of an edge-on disk whose major axis is perpendicular to the axis 
of the bipolar lobes. The MIDI spectrum and the visibilities of Mz 3 exhibit a clear signature of amorphous silicate, in contrast to 
the signatures of crystalline silicates detected in binary post-AGB systems, suggesting that the disk might be relatively young. We 
used radiative-transfer Monte Carlo simulations of a passive disk to constrain its geometrical and physical parameters. Its inclination 
(74°±3°) and position angle (5°+5°) are in accordance with the values derived from the study of the lobes. The inner radius is 9+lAU 
and the disk is relatively flat. The dust mass stored in the disk, estimated as 1 x 10"^ Mq, represents only a small fraction of the dust 
mass found in the lobes and might be a kind of relic of an essentially polar ejection process. 

Key words. Planetary Nebulae-Individual: Mz3 Techniques: interferometric; Techniques: high angular resolution;Stars: circumstel- 
lar matter; Stars: mass-loss 



1. Introduction 

■ It is commonly accepted that disks can be an essential ingredi- 
, ent to the shaping of planetary nebulae, but the spatial resolu- 
' tion of astronomical instruments is usually not sufficient for de- 
tecting and studying these disks. The geometry of the disk and 
, the mass stored are key parameters for constraining the mod- 
' els of nebula formation, and for tracing back the evolution of 
the central star CBahck & Frank 2002 ). Its inner edge, as seen 
from the star, could be thick and dense enough to collimate 
stellar winds into lobes, and knowing its geometry enables us 
to better understand the distribution of illuminated and shad- 
owed regions in the extended nebula. The VLT Interferometer 
is a powerful 'disk hunter' with its two instruments: AMBER 
operating in the near-IR dPetrov et al. 2007l l and MIDI in the 
mid-IR (Leinert et al. 2 003 ) providing a typical spatial resolu- 
tion of 2 and lOmas, respectively. It must be stressed that the 
number of visibility measurements usually recorded do not al- 
low any image reconstruction from the interferometric obser- 
vations because this would be too time-consuming. Therefore, 
the interpretation uses models that depend critically on the 
information already gathered on the sources by other tech- 
niques. As an example, a disk was detected by HST and MIDI 
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in the center of the complex multipolar nebula CPD-56°8032 
(IDe Marco et al. 20021 IChesneau et al. 2QQ6i . The disk exhibits 
a bright and extended (on an interferometric scale, ~70mas) in- 
ner rim whose PA and inclination do not seem directly related 
to extended structures in the HST images. Definitely, the disk 
detection based on three visibility measurements provided a lot 
of information, but also demonstrated the need for an interfer- 
ometric imaging campaign involving wide uv coverage for this 
complex object. 

In contrast, the Butterfly (M2-9) and the Ant (Mz 3, Hen 2- 
154) nebulae were chosen as perfectly suited targets to a snap- 
shot study with the VLTI, as they are among the most impres- 
sive bipolar PNs, highly collimated on large scales with tightly 
pinched waists as seen in Hubble Space Telescope (HST) im- 
ages (Balick & Frank 2002 and references therein; [Smith 20031 
Smith & Gehrz 2005. Smith et al. 2005J . For these targets, as- 
suming that the putative disk hidden in their center shares 
the symmetry of the extended nebula constitutes a good start- 
ing point for preparing and interpreting the VLTI observations. 
The nebula Mz 3 shows complex structures revealing different 
stages of mass-loss events that share the same axis of symme- 
try jSantander-Garcia et al. 20041 IGuerrero et al. 2004l l. Of im- 
portance is also the fact that M2-9 and Mz 3 are spectroscopic 
twins at visual and near-IR wavelengths. Each nucleus shows 
a rich optical spectrum with the various iron and helium lines, 
that are usually associated with a dense circumnuclear cocoon 
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which reprocesses most of the stellar light emitted in our line 
of sight that hides any close companion from direct view. Each 
object requires an ionization source with a temperature on the 
order of 30,000 K, and the dust mass detected in their nebulae 
suggests intermediate-mass progenitors (.Smith & Gehrz 2005I I. 
Published distances to each object vary widely, hovering around 
1-2 kpc. A detailed imaging study of the dust in the core of these 
nebulae has been presented recently by Smith (2003) and Smith 
& Gehrz (2005). The spatial resolution of the mid-IR images was 
1", allowing them to isolate the central engines from the sur- 
rounding nebulae. The observed mid-IR SED of the unresolved 
central sources is flatter than a single gray-body, implying that a 
range of dust temperatures and an additional component of hot 
dust are both required to fit the near-IR continuum spectrum. 
This suggests the existence of a circumstellar disk. 

The study of Smith & Gehrz (2005) motivated us to use the 
unique MIDI/VLTI capabilities for detecting the disks in the cen- 
ter of Mz 3 and M2-9. The present Letter is entirely devoted to 
the disk detected in the center of Mz 3, while M2-9 observations 
are described in another paper. The observations are presented 
in Section 2. In Section 3 we derive the physical parameters of 
the dusty environment by radiative-transfer models and discuss 
the results in Section 4. 

2. Observations 

Mz 3 was obs erved with MIDI dLeinert et al. 20031 
IRatzka et al. 2007l l the mid-infrared recombiner of the Very 
Large Telescope (VLT). The MIDI/VLTI interferometer op- 
erates like a classical Michelson interferometer to combine 
the mid-IR light (N band, 7.5-13.5jum) from two VLT Unit 
Telescopes (8.2m, UTs) or two Auxiliary Telescopes (1.8m, 
ATs). In our case, only the UTs were used for detecting 
correlated fluxes as low as about 6Jy from the 30Jy cores of 
Mz 3 and M2-9. We used a typical MIDI observing sequence, 
as described in Ratzka et al. (2007). MIDI can make single-dish 
acquisition images with a field of view of about 3" with a 
spatial resolution of about 0.25" at 8.7//m, and provides a 
flux-calibrated spectrum at small spectral resolution (R=30 
in this case) and several visibility spectra from the sources 
jChesneau et al. 2005al IChesneau et al. 2005bi) . i.e. spectrally 
dispersed information on the spatial extension of the source. 
The MIDI observations of Mz 3 were performed in May and 
June 2006 in the so-called SCIJPHOT mode, meaning that the 
photometry of each telescope is recorded simultaneously to the 
fringes. The errors, including the internal ones and the ones 
from the calibrator diameter uncertainty, range from 3% to 7%. 
The accuracy of the absolute flux calibration is about 10-13%. 
The log of the observation is given in Table [T] We used two 
different MIDI data reduction packages: MIA developed at the 
Max-Planck-Institut fiir Astronomic and EWS developed at the 
Leiden Observatory (MIA+EWSQ ver 1.5.1). 

We also observed Mz3 with the near-infrared AO instru- 
ment NACO using three broad-band filters centered respectively 
at 2.18//m (Ks), 3.8//m (L'), and 4.78//m (M'). The Ks and L'- 
band data were taken during the night of 18 August 2006 and 
the M'-band data were taken during the night of 29 May 2006. 
These observations were complemented with nearby PSF cali- 
brations using HD 155078 (F5V) in L' and HD156971 (FlIII) 
in M'. We used the S27 (Ks) and L27 (L',M') camera mode to 
obtain a field of view of 28"x28" and the pixel scale is 27. Imas 
per pixel. Jittering was used to remove the sky, and chopping was 

' Available at |http://www.strw.leidenuniv.ni/~nevec/MIDI/index.html| 
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also used for the M' observations. Data reduction was performed 
using ESO pipeline and self-developed IDL routines. 

Figure [T] shows the ISO and MIDI spectra. The ISO spec- 
tra are extensively discussed in Pottasch & Surendiranath (2005), 
and their respective apertures are shown in their Fig.l. We only 
show the SWS observations with the 14"x27" aperture. The flat 
MIDI spectrum contrasts to the steeply rising ISO flux toward 
the long wavelengths. Another obvious difference between the 
two spectra is that the 12.8/im [Nell] line is barely visible in the 
MIDI spectrum. 

The MIDI flux agrees within the 10% error bars with the 
fluxes reported with 4" apertures by Smith & Gehrz (2005, Table 
2) at 8.9 and 1 1.6yL(m, but the MIDI flux is significantly lower in 
the [Nell]. The M' band magnitude of the core agrees within er- 
ror bars with the ISO flux showing that the whole M' band flux 
originates in the core (mM'=L9+0.1). In contrast, the L' band 
flux of the core is much lower than the ISO curve (m^,- =4.4+0.1, 
aperture 0.4"), due to a strong reddening occurring at this wave- 
length and the increase in the light diffusion in the lobes. This 
magnitude value follows the trend in the fluxes measured with 
decreasing apertures in Cohen et al. (1978). 



3. Physical parameters of the disk 

The NACO and MIDI images show that the core is unresolved 
from the Ks band on a scale of about 60mas to the N band on a 
scale of about 250mas. The large amplitude of variations of the 
MIDI visibilities versus the baseline position angles can imme- 
diately be interpreted in the frame of a very flat structure, whose 
symmetry axis is roughly aligned with that of the bipolar lobes. 
In Table[T] and in Fig.|2l the projected baselines are organized 
from the ones providing the highest visibilities (i.e. smallest ex- 
tension) at PA close to 0° to the ones with lowest visibilities (i.e. 
highest extension) distributed around PA=90°. A rough estimate 
of the flux distribution using Gaussian ellipses provides sizes of 
lOmas X 20mas at 8yL(m, and 20mas x 40mas at 13;um. Another 
particularly interesting aspect of these visibility curves is their 
amorphous silicate signature, which is less obvious in the MIDI 
spectrum and hardly seen in the ISO spectrum. A higher opacity 
at a particular wavelength will cause a larger radius to be ob- 
served: the interferometer will detect this as a reduced visibility 
as compared to the continuum emission. This well-known opac- 
ity effect is particularly striking in the case of OH/IR stars (see 
IChesneau et al. 2005bl l. In the case of Mz 3, the similar shape 
of the visibility curves allows us to attribute it without ambi- 
guity to amorphous silicate. Mz 3 is a well-studied PN, and the 
a-priori information on its large scale geometry and dust content 
severely limits the range of spatial and physical parameters of 
the detected disk. The number and quality of the interferometric 
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Table 2. Model parameters and those from the model fitting. 



data were enough to conduct a direct interpretation via a radia- 
tive transfer model of a passive disk. 

3.1. A-priori knowledge on the disk of Mz3 

The large-scale spatial information originates essentially in 
the extended bipolar lobes (Santander-Garcia et al. 2004; 
IGuerrero et al. 2004l l. whose kinematic age is limited to less than 
ISOOyrs. The estimated inclination of the different structures 
ranges from 68° to 78°, and lGuerrero et al. 2004l mention a PA of 
10+2°for the lobes. The mid-IR study from lSmith & Gehrz 20051 
also provides some important information: the mass of warm 
dust inside the core is rather limited (minimum of 1 . 1 x IO^^Mq, 
optically thin environment and small carbon grain hypothesis), 
and the ratio between the mid-IR and visible luminosity suggests 
a modest opening angle (<24°). 

3.2. A passive disk model 

We made use of the 3D continuum radiative transfer code MC3D 
(Wolf 2003; see also Wolf et al. 1999). It is based on the Monte- 
Carlo method and solves the radiative transfer problem self- 
consistently. MC3D is designed to simulate dust temperatures in 
arbitrary dust/electron configurations and the resulting observ- 
ables: spectral energy distributions, wavelength-dependent im- 
ages, and polarization maps. We use a classical model of a dusty 
stratified disk (Shakura & Sunyaev 1973; Wood et al. 2002). The 
dust density follows a 2D law (radial and vertical dimensions): 



p(r, z) = po 



exp 



z 

W) 



where r is the radial distance in the midplane of the disk, a the 
density parameter in the midplane, 7?* the stellar radius and the 

disk scale height h{r) is given by h{r) - /lo(^) , where is 
the scale height at a given radial distance from the star and fi is 
the vertical density parameter. 



We assume the standard interstellar grain size distribution 
(Mathis et al. 1977): a a"-''^ where a is the dust grain ra- 
dius. Grain radii extend from 0.05 to Ifim and we consider the 
dust grains to be homogeneous spheres. The MIDI data are not 
suited to constrain the dust size distribution efficiently, so that the 
dust grain sizes are assumed and kept fixed in the fitting process. 
Similarly, the outer radius of the structure is not well-defined by 
these N band observations, so it is arbitrarily limited to 500 AU 
(0.35" at 1.4 kpc). The outputs from the code are a large aper- 
ture SED, 15 images (2mas per pixel) in the wavelength range 
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Fig. 1. MIDI spectrum (thick line) compared with the ISO spec- 
trum (upper curve) and the TIMMI2 photometric measurements 
of Smith & Gehrz (2005) with an aperture of 4" (triangles). The 
NACO L' and M' bands are also included (stars, aperture 0.4"). 
The dashed line is the best silicate model with the full aperture 
(1"). In the N band, the 1" beam and the MIDI beam flux of the 
model differ by only a few percent. 

8-13.5 jum, and by applying 2D Gaussian apertures to the im- 
ages (FWHM from 210mas to 350mas), we get the MIDI spec- 
trum and the spectrally dispersed visibilities (taking the limited 
spectral resolution into account). To generate the visibility sig- 
nal, the images are first rotated by the PA of the object in the 
sky and then collapsed as ID flux distributions in the direction 
of the baselines. The visibility for each wavelength is the value 
of the Fourier transform of these ID vectors, normalized to the 
zero frequency value. 

The distance, luminosity, and temperature of the central star 
being uncertain by a large factor, we fixed them to plausible val- 
ues close to the ones in Sm ith 20031 trying to find an adequate 
fit of the SED and visibilities with them. It happened that such 
a good fit was found rapidly, and we did not try to explore the 
parameter space for the central source further, because we think 
that the current data are not sensitive enough to distinguish one 
model from another. With such a large inclination, the fi and h 
parameters are well-fitted using the balance of visibilities be- 
tween baselines perpendicular and parallel to the structure. 

4. Results and discussion 

The MIDI observations provide evidence of a flat, nearly edge- 
on disk primarily composed of amorphous silicate whose pa- 
rameters are described in Table 2. The inner rim is at sublima- 
tion temperature and the disk is optically thick {t^ ~ 3.5). The 
dust chemistry agrees with |Bernard-Salas & Tielens (2005) ] and 
[Cohen & Barlow (2005)| who found evidence in the ISO spec- 
tra of weak signatures of amorphous silicate and no obvious 
signature of carbon dust in the form of PAHs, in accordance 
with the CIO ratio of 0.83 determined by [Zhang & Liu 2002 
in the gas phase. The match of the SED and the visibil- 
ity curves is impressively good showing that such a simple 
model of a passive disk can account for these observations. 
That the silicate is in the form of amorphous grains is impor- 
tant, since the long-lived disks found in binary post-AGBs ex- 
hibit strong signatures of crystalline silicate dDeroo et al. 2007! 
De roo & Van Winckel 2 007). This suggests that the disk may be 
relatively young, in line with the youth of the PN. The dust mass 
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Fig. 2. MIDI visibilities with their error bars, compared with the best silicate model (dashed-dotted lines, see Table |2]i. The visibility 
curves are Mz3-1 to Mz3-6 from top to bottom and left to right, with the;t'^ values of the fits of 2.16, 0.30, 0.95, 1.50, 0.88, 4.31. 
The reduced for the full data set is 1.68. The more extended [NeII]12.8jL/m line is not taken into account in the model. 



stored in the disk, estimated to be 1 x 10"^ Mq is two orders of References 
magnitude below the mass inferred in the lobes (2.6 x IO^^Mq 
from ISmith & Gehrz 20051 1. The low mass of the disk seems 
insufficient for pinching the waist of the bipolar nebula hy- 
drodynamically, pointing instead to an essentially polar ejec- 
tion process. The small disk opening angle and the limited an- 
gle sustained by the lobes also favor such a scenario. The de- 
tection of a bright X-ray core and a possible jet, reported by 
IKastner et al. 20031 is probably a good indication that a compact 
accretion disk is hidden from direct view by the more extended 
dusty disk presented here. Using near-IR color-color diagram 
and the similarity between the locus of symbiotic stars and Mz 3, 
|Schmeja & Kimeswenger (2001) 1 suggested that a Mira should 
be hidden in the center of the system. Our preferred interpre- 
tation is that this kind of diagram is useful for probing dusty 
disks, but is not convincing evidence of a Mira star. Only the 
temporal variability of the near-IR flux would represent such ev- 
idence but it is currently lacking for Mz 3. The constraint pro- 
vided by our observations is that the orbit of the system must be 
well within the inner radius of the disk, at a distance of at most 
1-2 AU from the primary, which is typical of S-type symbiotics 
(Mikolajewska 2003 1 but much too small for D-type symbiotic 
Miras. Therefore a Mira seems improbable, but a less luminous 
cool giant companion might remain hidden by the disk. The con- 
straints from the SED are tight, limiting the cool companion lu- 
minosity to less than IOO-I5OL0 (ISmith 20031) . This object could 
be a young link to the binary post-AGBs ( de Ruyter et al. 2006[ l 
that probably experienced non conservative mass loss. The cur- 
rent estimate of the mass stored in the Mz 3 disk is in line with 
the mass of long-lived (crystalline) disks found in many of the 
binary post-AGBs. 
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